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ABSTRACT 

Context. Identifying the main processes of galaxy assembly at high redshifts is still a major issue to understand galaxy formation and 
evolution at early epochs in the history of the Universe. 

Aims. This work aims to provide a first insight into the dynamics and mass assembly of galaxies at redshifts 1.2 < z < 1.6, the early 
epoch just before the sharp decrease of the cosmic star formation rate. 

Methods. We use the near-infrared integral field spectrograph SINFONI on the ESO-VUT under 0.65" seeing to obtain spatially 
resolved spectroscopy on nine emission line galaxies with 1.2 < z < 1.6 from the VIMOS VUT Deep Survey. We derive the velocity 
fields and velocity dispersions on kpc scales using the Ho- emission line. 

Results. Out of the nine star-forming galaxies, we find that galaxies distribute in three groups: two galaxies can be well reproduced 
by a rotating disk, three systems can be classified as major mergers and four galaxies show disturbed dynamics and high velocity 
dispersion. We argue that there is evidence for hierarchical mass assembly from major merger, with most massive galaxies with 
M > 10" Mq subject to at least one major merger over a 3 Gyr period as well as for continuous accretion feeding strong star 
formation. 

Conclusions. These results point towards a galaxy formation and assembly scenario which involves several processes, possibly acting 
in parallel, with major mergers and continuous gas accretion playing a major role. Well controlled samples representative of the bulk 
of the galaxy population at this key cosmic time are necessary to make further progress. 
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0^ 1 . Introduction narl o, whereby the mo st massive galaxies assemble their mass 

O . . . . , first ^^EmMMBi- 
In the current paradigm, galaxies are formed in dark matter 

. !^ haloes which grow along cosmic time via hierarchical assem- However, much remains to be done to understand how galax- 
^ . bly of smaller units. While in this picture halo-halo merging is ies of different types have been built along cosmic time. Several 
^ I the main physical process driving the assembly of dark matter key questions are still subject to considerable debate: (i) How 
5^ ' haloes, it is unclear how this merging process directly affects the were galaxies assembled in details? (ii) When and at what rate 
build-up of galaxies, and how other physical processes play a did galaxies of different masses form? (iii) What is the connec- 
role. tion between bulge and disk formation? (iv) What is the link 
Considerable efforts have been invested to understand galaxy between the various high-z galaxy populations, their origin and 
evolution in the past two decades. We now have a global, but their subsequent fates? (v) What is the connection between star 
as yet incomplete, picture of the evolution of the parameters de- formation and the AGN phenomenon ? 
scribing the main galaxy population, with e.g. the gl obal star for- 
mation rate reaching a peak at re dshifts 1 to 2.5 (e.g. lTresse etalJ Several physical processes have been identified to be con- 
l2007tlHopkins & Beacoml2006l) . A clear result now emerging is tributing to the evolution of galaxies, but their relative contri- 
that early type galaxies have experienced a major growth in stel- butions, main epoch of action, and associated timescales, re- 
lar mass density from a redshift z ~ 2 to z ~ 1, while late-type main poorly constrained. Merging is identified as an important 
galaxies hav e been growing their stellar mass more sl owly in this contributor through direc t evidence of on-g oin g events, o rjndi- 
period re.g. lArnouts et alJ 120071: iBundv et alJl2005l) . It is also rect fossil remnants (e.g.|Conselice et aU2003HLotz et al.|2008|), 



claimed that galaxy evolution proceeds via a downsizing see- and the merger rate and its evolution is now demonstrated to 

be strongly dependent on the luminosity o r stellar mass of the 

* Based on observations collected at the European Southern galaxies involved (e.g. |de Ravel et aljUOOl). In addition, the ef- 

Observatory (ESO) Very Large Telescope, Paranal, Chile, as part of the feet of the competing processes of cooling, angular momentum 

Programs 75.A-0318 and 78.A-0177 exchange, feedback from star format ion and AGNs, or gas ac- 

** E-mail: benoit.epinat@ast.obs-mip.fr cretion, are to be further investigated (Somerville et alfcOOSi) . 
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Until recently, the favoured picture was that the star for- 
mation and stellar mass assembly of massive galaxies was a 
direct co nsequence of major m erging events at early cosmic 
epochs ( Glazebrook et al .1120041) . However, recent observations 
of individual galaxies have revealed that large and massive 
disks with strong star formation seem to be already in place 
at redshifts z ~ 2, without apparent signs for major merging 
events ( Forster Schreiber et al. 2006; Wright etal. 2007, 2009; 
iLaw et al . 2007, 2009; Genzel et al. 2006, 2008). This has led to 
suggest that the majority of star forming galaxies is fed by gas 
via cold flows along streams of the cosmic web (in cluding mi- 
nor m erger), continuously fueling the star formation (iDekel et al.l 
12009) . This scenario is further supported by the tight correla- 
tion between stellar mass and star forrnation rate in high redshift 
star forming galax ies (e.g.' Paddi et alj2007tlNoeske et al.l2()07t 
lElbaz et al.ll2007l) . Dynamical processes internal to galaxies can 
then drive secular evolution of disks and the formation of bulges 
and spheroids. 

The detailed knowledge of the kinematics of a galaxy on 
kiloparsec scales is needed to identify its dynamical state (disk 
in rotation, spheroid, major marger event or more complex kine- 
matics). The integral field technique enables one to compute ac- 
curate total dynamical masses, to trace the spatial distribution of 
stars and gas, as well as to evaluate the contribution of stellar 
populations including stellar initial mass function (IMF), or gas 
metallicity. The velocity field of spiral/disk galaxies, can be used 
to put important constraints on total masses and hence on dark 
matter halo masses, thought to be an important driver of galaxy 
evolution. 

At low redshifts, the knowledge of 2D velocity fields has 
proved to be a powerful kin ematical tool to i nvestigate the 
prope r ties of nearby galaxies jOstlin etaTI 120011; IVeilleux et al.l 



2001'; 'Swinbanket al.' '20 031; iMendes de Oliveira et all 120031; 



Garrido et al. 2004; Epinat e t al .Il2008allbh . The velocity fields of 
galaxies at intermediate redshift (0.4 < z < 0.75 ) have been 
inves t igated using FLAMES/GIRAFFE at VLT dPlores et alj 



mves t igateq usmg i^LAMi^,t)/uiKA^i^i^, at vi^i Qi^iores et all 
2006t liPuech et alJ 120061: IVang et al.ll200i iNeichel et al.ll2008t 



Puech et al. ^2008). In their sample of 65 galaxies. lYang et akl 
(2008) found about 32% of relaxed rotating disks. These rotat- 
ing disks produce a Tully-Fisher relation tTullv & Fisheij[T977h 
which has apparently not evolved in slope and scatter since 
z=0.6. They however detect an evolution of the K-band Tully- 
Fisher relation zero point that, if interpreted as an evolution of 
the K-band luminosity, would imply a brightening of 0.66±0.14 
mag between z ~ 0.6 and z = 0. They suggest that the large scat- 
ter found in previously reported Tully-Fisher relations at moder- 
ate redshifts are produced by the numerous (65%) galaxies with 
perturbed or complex kinematics. 

At higher redshifts the SINS survey has pioneered the de- 
tailed observations of galaxies at z ~ 2 (jp orster Schreib er et al 



2006 1: iGenzel et alJl2006ll2008tlShapiro et al...2008i;iCresci et al 



20091) . demonstrating that near infrared integral field spec- 
troscopy is the best tool to securely measure the kinematical 
properties at high redshifts z > 1. These authors have found 
that galaxies have an extremely large velocity dispersion (cr ~ 
80 km 5"') as compared to their rotational velocity. They fa- 
vor the hypothesis of early buildup of central disks and bulges 
by secular evolution in gas-rich disks. Fast turbulent speeds in 
the gaseous component imply the formation of massive clumps. 
Gas-rich primordial disks may evolve through a clumpy phase 
into bright early-type disk galaxies with a massive exponen- 
tial disk, a classical bulge and possibl y a central black hole 
(Noguchi 1999; Immelietal. 2004a ibt Elmegreen et al] 120051: 
iBournaud et al.i i2007t, ,2008 ). These massive clumps are sug- 



gested from NACO high resolution deep imaging and Ha maps 
from adaptive optics SINFONI observation in two galaxies of 
the SINS sample (Genzel et al. 2008). 

The number of galaxies observed at or around the peak in 
star-formation at z ~ 1 - 2.5 is still relatively small, and detailed 
observations of a large volume complete sample of galaxies se- 
lected from an homogeneous set of selection criteria are needed 
to clarify the respective contributions of merging and continuous 
gas feed by accretion. 

Assembling large samples of galaxies is not an easy task as 
we need to extract from large spectroscopic surveys, selected 
with well controlled criteria, those galaxies with accurate spec- 
troscopic redshifts for which the Hg line emission falls outside 
bright OH sky lines. The VVDS dLe Fevre et alj|2005h survey 
has been designed to map the evolution of galaxies, large scale 
structures and AGNs from the spectroscopic redshift measure- 
ment of tens of thousands of objects. This dedicated program 
has successfully crossed the "redshift desert" 1 < z < 3, pro- 
viding for the first time a complete magnitude-limited sample of 
~ 800 galaxy redshifts in the "Ultra-deep" (Iab ^ 24.75) sample 
(Le Fevre et al., in preparat ion), ~ 11 000 galaxy redshifts in the 
"deep" (Iab < 24) sample dLe Fevre et alj|2005l). and ~ 40000 



in the "wide" (Iab < 22.5) sample dGarilh et al.ll2068h . 

We are using VVDS to select targets for a large observing 
program at the ESO-VLT aimed at probing the mass assembly 
and metallicity evolution of a representative sample of galax- 
ies at 1 < z < 2, a crucial epoch corresponding to the peak of 
cosmic star formation activity. The main goal of the MASSI\Q 
(Mass Assembly Survey with SINFONI in VVDS) project is to 
obtain a detailed description of the mix of dynamical types at 
this epoch and to follow the evolution of fundamental scaling 
relations, such as the Tully-Fisher or the mass-metallicity rela- 
tions, and therefore constrain galaxy evolution scenarios. 

In this paper we present the first results obtained by MASSIV 
focusing on the kinematics of nine galaxies with 1 .2 < z < 1 .6 
observed during the MASSIV pilot program. A companion pa- 
per (Ouevrelet al. 2009) is devoted to the mass-metallicity re- 
lation of galaxies at these redshifts using the same data. In (§2) 
we describe our observations and data reductions. In (§3) we 
present the morphological, physical, kinematical and dynami- 
cal measurements of the galaxies of our sample. In (§4 ) we 
classify the galaxies exploiting the full kinematical information 
and discuss the Tully-Fisher relation. Our discussion and conclu- 
sions are provided in (§5). Appendix lAl contains detailed indi- 
vidual comments for each galaxy. We assume a cosmology with 
0,„ = 0.3, Q.A = 0.7 and //q = 70 km s"' Mpc"' throughout. 



2. Data and Observations 

2.1. Sample selection 

We have used the VIMOS VLT Deep Survey sample to select 
galaxies with known spectroscopic redshifts across the peak of 
star formation activity. The VVDS is a complete magnitude se- 
lected sample avoiding the biases linked to a priori color selec- 
tion techniques. This sample offers the advantage of combining a 
large sample with a robust selection function and secure spectro- 
scopic redshifts. The latter are necessary to engage into long sin- 
gle objects integrat ions with SINFONI dEisenhauer et alj|2003l; 
iBonnet et al.l 120041) being sure to observe the Ha line outside 
bright OH night-sky emission lines. 



www. ast . obs-mip . fr/massiv/ 
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Table 1. Table of observations 
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VVDS ID a (5 z /ab Field t„p Seeing Run ID 

J2000 J2000 mag hours 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



020116027 


02 


25:51.133 


-04:45:04.48 


1.5259 


22.875 


VVDS-02h 


1.67 


0.65 


075.A-03 18(A) 


020182331 


02 


26:44.260 


-04:35:51.89 


1.2286 


22.729 


VVDS-02h 


3 


0.72 


078.A-0177(A) 


020147106 


02 


26:45.386 


-04:40:47.39 


1.5174 


22.502 


VVDS-02h 


2 


0.74 


075 .A-03 18(A) 


020261328 


02 


27:11.049 


-04:25:31.60 


1.5291 


23.897 


VVDS-02h 


1 


0.61 


075.A-0318(A) 


220596913 


22 


14:29.184 


-1-00:22:18.89 


1.2667 


21.841 


VVDS-22h 


1.75 


0.47 


075 .A-03 18(A) 


220584167 


22 


15:23.038 


+00:18:47.01 


1.4637 


22.036 


VVDS-22h 


1.75 


0.74 


075 .A-03 18(A) 


220544103 


22 


15:25.708 


+00:06:39.53 


1.3970 


22.469 


VVDS-22h 


1 


0.69 


075 .A-03 18(A) 


220015726 


22 


15:42.455 


+00:29:03.59 


1.3091 


22.473 


VVDS-22h 


2 


0.58 


075 .A-03 18(A) 


220014252 


22 


17:45.690 


+00:28:39.47 


1.3097 


22.101 


VVDS-22h 


2 


0.70 


075 .A-03 18(A) 



(1) Source VVDS identification number, (2) and (3) Coordinates, (4) VVDS spectroscopic redshift, (5) AB magnitude in I-band, (6) VVDS-22h 
wide field (17.5 < /ab < 22.5) and VVDS-02h deep field (17.5 < /ab < 24.0), (7) Exposure time, (8) Median seeing of SINFONI observations 

measured on PSF stars, (9) ESO program. 
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Fig. 1. Selection of late-type 1 < z < 1.54 star-forming galaxies 
with secure redshift and measured [O ii] /13727A emission line 
in the VVDS-22h (/ab < 22.5, small triangles) and VVDS-02h 
(/ab ^ 24, small squares) fields for SINFONI follow-up obser- 
vations. Among these galaxies we selected those showing the 
strongest [On] /13727A emission line (filled symbols: EW > 50 
A and flux > 5 x 10"'^ erg s^'cm"^, measured on VIMOS spec- 
tra) for Ha to be easily detected in the near-IR. The nine galaxies 
studied in this paper are indicated as large blue triangles (VVDS- 
22h field) and large red squares (VVDS-02h field). 



From the existing VVDS dataset, we have available today 
a unique sample of more than 1500 galaxies in the redshift do- 
main 1 < z < 2 with accurate (to the fourth di git) and secure 
spec troscopic redshifts (confidence level > 90%: lLe Fevre et al] 
12005 ). This sample being purely /-band limited, it contains both 
star-forming and passive galaxies distributed over a wide range 
of stellar masses, enabling to easily define volume limited sub- 
samples. For the MASSIV program, we have defined a sample 



of ~ 140 VVDS star-forming galaxies at 1.0 < z < 1.8 suit- 
able for SINFONI observations. In most of the cases, the galax- 
ies are selected on the basis of their measured intensity of [O ii] 
A3727 emission line in the VIMOS spectrum or, for a few cases, 
on their observed photometric U BVRIK spectral energy distri- 
bution which is typical of star-forming galaxies. The star for- 
mation criteria ensure that the brightest rest-frame optical emis- 
sion lines, mainly Ha, [Oiii] ^^4959,5007, [Nii] /16584 used 
to probe kinematics and chemical abundances, will be observed 
with SINFONI in the NIR bands. 

For the pilot observations presented in this paper, the selec- 
tion has been made on [O ii] A3121 emission line only thus re- 
stricting the redshift range to 1 < z < 1 .54. We have selected 
galaxies showing the strongest [On] A3727 emission line (EW 
> 50 A and flux > 5 x 10"'^ erg s^'cm"^) as measured on 
VIMOS spectra (see Figure [T]) for Ha to be easily detected in 
the near-IR. Among these candidates, we have further restricted 
the sample taking into account one important observational con- 
straint: the observed wavelength of Ha line had to fall at least 
lOA away from strong OH night-sky lines, to avoid heavy con- 
tamination of the galaxy spectrum by sky subtraction residuals. 
Among the remaining list of promising candidates, we have ran- 
domly picked up twelve galaxies to carry out these first pilot 
observations. These criteria for selecting late-type star-forming 
galaxies have been shown to be very efficient. For the first ob- 
serving runs, our success rate of selection has been around 85%: 
9 galaxies over 12 observed show strong rest-frame optical emis- 
sion lines in SINFONI datacubes, the 3 remaining show a signal 
too faint to derive reliable velocity maps. Among the nine VVDS 
star-forming galaxies, five have been selected in the VVDS- 
22h (17.5 < /ab < 22.5) field and four in the VVDS-02h 
(17.5 < /ab ^ 24.0) field. These targets span a redshift range 
between z ~ 1.22 andz ~ 1.53 (see Table[T]). 

2.2. SINFONI observations and data reduction 

The NIR spectroscopic observations were obtained with the 
3D-spectrograph SINFONI at ESO-VLT during two four-nights 
runs, on September 5-8, 2005 (ESO run 75.A-0318) and on 
November 12-15, 2006 (ESO run 78.A-0177). SINFONI was 
used in its seeing-limited mode, with the 0.125"x0.25" pixel 
scale leading to a field-of-view of 8"x8", and the H grism pro- 
viding a spectral resolution R ~ 4000. Conditions were not pho- 
tometric and the mean seeing measured on PSF stars was around 
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Table 2. Fluxes and luminosities of VVDS galaxies 



VVDS ID Vab L,5oo FCHa) L(HQr) F([On]) L([Oii]) 

mag 10^^* ergs s"' Hz"' 10"'^ ergs s"' cm"^ 10^"*' ergs s"' 10"" ergs s"' cm"^ lO^*' ergs s"' 
(1) (2) (3) (4) (5) (6) (7) 



020116027 


24.0 


5.1 d 


:0.2 


6± 


1 


9± 


2 


6d: 


2 


9± 


2 


020182331 


24.1 


3.7 d 


:0.1 


55 ± 


14 


47 ± 


13 


20: 


t 2 


17 = 


t 2 


020147106 


23.4 


12.9 


±0.1 


22 : 


t6 


32: 


t9 


9d 


2 


12 = 


t3 


020261328 


24.2 


4.4 d 


:0.2 


5 ± 


1 


7 ± 


2 


10: 


t 2 


15 = 


t3 


220596913 


23.1 


8.7 d 


I 0.4 


33 : 


t6 


32: 


t6 


19: 


t 2 


17 = 


t 2 


220584167 


22.9 


12.5 


±0.6 


48 ± 


10 


64 ± 


14 


18: 


t 2 


24 = 


t3 


220544103 


23.3 


9.4 d 


:0.5 


161 : 


t29 


192: 


t 35 


25 : 


t3 


30 = 


t4 


220015726 


24.0 


3.1 d 


I 0.4 


26: 


t3 


26: 


t 3 


9d 


2 


10 = 


t 2 


220014252 


22.7 


11.2 


±0.3 


212 : 


t27 


216: 


t28 


38: 


t 2 


38 = 


t 2 



(1) Source VVDS identification number, (2) AB magnitude in (/-band, (3) UV luminosity, (4) SINFONI Ha flux, (5) Ho- luminosity, (6) VIMOS 

[On] /13727 flux, (7) [On] ^3727 luminosity. 



Table 3. SFRs and extinctions of VVDS galaxies 



VVDS ID 


SFRt/v 


SFRhii 


SFR[o„] 


SFRsED 


E(B - V)uVIHa 


E{B - V)SED 




Mo yr-' 


Mo yr-' 


Mo yr-' 


Mo yr-' 


mag 


mag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


020116027 


7.2 ± 0.3 


7 ± 2 


13 + 4 


46 


0.0 


0.2 


020182331 


5.2 ± 0.2 


38 ± 10 


24 + 7 


67 % 


0.2 


0.3 


020147106 


18.1 ± 0.2 


25 + 7 


17 + 6 


55 *l 


0.0 


0.1 


020261328 


6.2 ± 0.2 


6±2 


21 + 7 


8+f 


0.0 


0.1 


220596913 


12.1 ± 0.6 


24+5 


24 + 7 


77 +109 
-29 


0.1 


0.1 


220584167 


17.5 ± 0.8 


51 ± 11 


34+ 10 


158 „2 
46 


0.1 


0.3 


220544103 


13.2 ±0.7 


152 + 28 


42 + 13 


0.3 


0.3 


220015726 


4.3 ± 0.5 


21 + 3 


13 + 4 


86 +1 
84 


0.2 


0.2 


220014252 


15.7 ± 0.4 


171 + 22 


54+ 16 


0.3 


0.3 



(1) Source VVDS identification number, SFRs deduced respectively from (2) UV flux, (3) Ha flux, (4) [On] ^3727 flux and (5) SEDs, 
Reddening sufi^ered by the ionised gas respectively computed (6) from the UV/Ha ratio and (7) from SED modeling. 



0.65" (see details in Table [Tli. The PSF stars were observed by 
SINFONI in //-band for each pointing on any target. Each target 
was acquired through a blind offset from a nearby bright star. 
Each observation was obtained by nodding the position of the 
galaxy within the 8" x 8" SINFONI field-of-view, generally 
by locating the source in two opposite corners. This observa- 
tional procedure allows background subtraction by using frames 
contiguous in time, but with the galaxy in different locations. 
Moreover, the target was never located exactly at the same po- 
sition on the detector, a minimal dithering of 0.3" was required 
in order to minimize instrumental artifacts when the individual 
observations are aligned and combined together. The total on- 
source integration times are listed in Table [T] and range between 
1 and 3 hours. As observations were performed in visitor mode, 
the integration time has been adjusted on the fly depending on 
the signal to noise achieved for the Ha emission line. 

Data reduction has bee n performed w ith the ESO-SINFONI 
pipeli ne (version 1.7.1, see lAbuter et al.ll2006: Modi gliani et al.i 
l2007h . The pipeline subtracts the sky background from the tem- 
porally contiguous frames, flat-fields images, spectrally cali- 
brates each individual observation and then reconstructs the dat- 
acube. Individual cubes were aligned in the spatial direction by 
relying on the telescope offsets from a nearby bright star and 
then combined together. A flux calibration is mandatory in order 
to derive absolute parameters (e.g. star formation rate) from the 
flux measured in emission lines. To this end, each science ob- 
servation has been immediately followed by the observation of a 
telluric standard star 



3. Measurements 

3.1. Stellar masses, star formation rates and extinction 

Stellar mass M^tar and star formation rate (SFR) estimates for 
our sample were obtained using the GO SSIP spectral energy 
distribution (SED) modeling software ( Franzetti et alJ l2008h . 
We used as input for the SED fitting the multi-band photo- 
metric observations available in VVDS fields, including BVRI 
data from the CFHT, UBVRZs data from the CFHT Legacy 
Survey, J- and /T-bands data from SOFI at the NTT and from 
the UKIDSS survey, and the VVDS-Deep spectra. The photo- 
metric and spectroscopic data were fitted with a grid of stellar 
population models, gen erated u sing the BC03 population syn- 
thesis code ( Bruzual & Charlol l2003i). assuming a set of "de- 
layed" star formation h istories (see iGavazzi et al.ll2002l for de- 
tails), a' SalpeteJ (Il955h IMF with lower and upper mass cutoffs 
of respectively 0.1 and 100 Mq, a metallicity ranging from 0.02 
and 2.5 solar meta llicity and a Large Magellanic Cloud redden- 
ing law ftPeil 19921) with an extinction E(B - V) ranging from to 
0.3. The parameters for the best-fitting model for each galaxy are 
taken as the best fitting values for both the galaxy stellar mass 
and SFR. On top of the best fitting values, GOSSIP computes 
also the Probability Distrib ution Function (PDF), following the 
method described in .Walcher et al.l (l2008h . The median of the 
PDF and its confidence regions are then used to derive a robust 
estimate of value and error for the parameter that is to be deter- 
mined. Stellar masses for the galaxies of this paper are reported 
in Table |5] The best SED fit model for each galaxy is presented 
in Figure |2l 
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Fig. 2. Best SED fit model BC03 (continuous line) and photometric points overplotted for the nine galaxies of the sample. 



In addition to the determination from the SED fitting 
(SFR^Eo), SFRs for the nine star-forming galaxies in our sam- 
ple can be deduced from the Ha (SFRetr) flux, from the [On] 
/13727 flux (SFR[Oi,]) and from the rest-frame UV continuum 
emission (SFRfyv). As these four estimators are affected differ- 
ently by dust and star formation history, our results can in prin- 
ciple tell us about the extinction and stellar populations of the 
galaxies. 

The Ha recombination line is a direct probe of the young, 
massive stellar population and therefore provides a nearly in- 
stantaneous (i.e. averaged on the last ten million years) measure 
of the SFR. Moreover, the Ha line is less aff'ected by dust extinc- 
tion com pared to e.g . UV e stimators. We have calculated SFRhc 
following lkennicutti ( Il998l) : 



SFRh. -.g,, 10-42 L(Ha) 
[Mayr-i] ' [ergs-'] 



(1) 



where the Ha luminosity L(Ha) ha s been derived from the Ha 
flux measured on SINFONI data by lOueyrel et"an (l2009l) . 

The [O ii] /13727 emission line luminosity L([0 ii]) is not di- 
rectly coupled to the ionising luminosity and its excitation is 
sensitive to abundance and the ionization state of the gas. It can 
however be empirically used as a quantitative SFR tracer. We 



have computed SFR[o„] following the calibration provided by 
Kennicutt (1998): 



SFR, 



[On] 



= (1.4 + 0.4) X 10 



-4iL([On]) 



(2) 



[Moyr-i] ^' ■ [ergs- 

where L([0 n]) has been derived from the [O n] /13727 flux mea- 
sured on VIMOS spectra. 

Ultraviolet-derived SFRs were calculated from the broad- 
band optical photometry. Given the redshifts of our targets, op- 
tical photometry provides indeed an information on the amount 
of rest-frame UV flux. The UV continuum around 1500 A is 
roughly flat if the flux is expressed in frequency units. The 
mean luminosity inside the optical broad-band filter which cor- 
responds to the rest-frame 1500 A continuum Lisoo, at the given 
redshift z of the source, is thus a very good approximation of the 
level of this continuum. The mean luminosity in frequency units 
is calculated following this equation: 

/ n. ^2 

Li500 -"-0.4 „ o^oi „ ^_V[ 



10 



[erg s 'Hz ' j 

Dl beeing the luminous distance. 
We have then deduced SFRj/y 

Ll500 



X 3631 x47r' 



- [cm] ; 
1 -HZ 



(3) 



=1.4xl0-2« 



[Mo yr-i] 



[erg s Hz ] 



followingl^rmicutd ( Il998l) : 

(4) 
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The three equations[T]|2]and|4]assume a Salpeter IMF with lower 
and upper mass cutoffs of 0.1 and 100 M© and do not include 
any effects of attenuation by dust. 

Table |2] gives the AB magnitude in f/-band, UV flux. Ho- 
and [On] /13727 fluxes and luminosities used to compute the 
four SFRs stored in Table |3] together with the reddening values 
deduced from SED fitting and from the UY/Ha SFR ratio. 

The UV SFR may be much more strongly affected by ex- 
tinction than the Ha one due to a shorter wavelength emission. 
Its comparison with the Ha SFR may thus provide an estimate of 
the amount of dust. Assuming that all of the ionizing photons are 
reprocessed into lines, therefore SFRhb = SFRf/y, we thus com- 
puted the interstellar gas extinction {E(B - V)) using the intrinsic 
starburst flux density fi(A) via: 

/;(^) = /„(i)io°-4^<^-^)*'w, (5) 

where the obscuration curve for the stellar continuum, k^{A), is 
given by Calzetti et al. (2000) for star-forming galaxies: 

k'(A) = 1.17 X (-2.156-H 1.509//1- 0.198// -hO.01 1/^^) + 1.78 



for O.llfjm < A < 0.63fjm 



(6) 



The amount of interstellar extinction estimated in our sam- 
ple (mean value E(B - V) ~ 0.15) is typical of star-forming 
galaxies in the local and intermediate-redshift universe and is in 
very good agreement with the reddening values derived at simi- 
lar redshifts in UV-selected galaxy samples by Erb et al. (2006). 
Figure |3] shows a good agreement between extinctions derived 
from SED fitting and from the UV/Hc SFR ratio. However, the 
SFR-based extinctions are systematically underestimated com- 
pared to the SED ones. This may be due to the uncertainties in 
the absolute flux calibration both in the rest-frame UV and even 
more important in the SINFONI data. 

In Figure |4] (left), we compare the SFR derived from 02 and 
Ha emission lines. The agreement is good taking into account 
that neither Ha nor 02 flux has been corrected for reddening. 
Indeed, Ha based SFR are systematically higher than 02 ones 
because 02 is more affected by extinction than Ha, except for 
the two galaxies with almost no extinction (see TableO. The fact 
that the two galaxies with the highest extinction (E(B-V) ~ 0.3) 
show the largest deviation is supporting this interpretation. 

The SFR deduced from the SED fitting (SFRs^b) provides 
an internal correction for dust extinction but is averaged on 
a timescale ten times longer than the SFR derived from Ha. 
The comparison between these two SFRs may provide informa- 
tion on the ratio between recent and older star formation activ- 
ity, and thus can give an hint on the presence of recent bursts 
(see Figure m right). Despite the large uncertainties inherent to 
the SED-based SFR determination, a good agreement between 
the two SFR estimators is seen for the three galaxies classified 
as perturbed rotators (see section 14.1b with a Ha-based SFR 
value lower than 50 Mq yr"'. The two galaxies classified as 
rotating disks and the two merging systems VVDS0201 16027 
and VVDS220596913 show an instantaneous Ha-based SFR 
(timescale ~ 10 Myr) which is lower (by a factor 2-4) than the 
SED-based SFR integrated over a longer timescale (a few hun- 
dreds Myr). On the contrary, the two galaxies with the highest 
values of Ha-based SFR (> 150 Mq yr"', VVDS220544103 and 
VVDS220014252) are typical starburst galaxies with a high ra- 
tio between the instantaneous and integrated SFRs. 

3.2. Morphologies 

For each galaxy of our sample we have estimated the inclination 
from the /-band CFHT images using the GALFIT software 
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Fig. 3. Comparison of extinctions computed from SED fitting (x- 
axis) and from the comparison of SFRs deduced from UV and 
Ha luminosities (y-axis). The dashed line indicates y = x. 



(lPengetal.1120021) . We have used CFHT le gacy survey images 
with best seeing (<0.65") for VVDS-02h field galaxies (Figures 
O to [AS and CFH12K/CFHT ( McCracken et all |2003|) 



images for VVDS-22h field galaxies (Figures IA.5I to |A.9l l. 
GALFIT convolves a point spread function (PSF) with a model 
galaxy image based on the initial parameter estimates fitting a 
ISersid d 19681) profile. Finally, GALFIT converges into a final set 
of parameters such as the Sersic index, the center, the position 
angle, and the axial ratio. Residual maps were used to optimize 
the results. Since the morphological inclination is used as input 
for kinematical modeling, we have imposed the morphological 
position angle of the major axis to be equal to the kinematical 
one (determined in section 13.3.2b . When the morphological 
position angle is a free parameter, the agreement with the kine- 
matical one is better than 25° except for galaxies which have a 
morphology compatible with a galaxy seen face-o n within the 
uncert ainties, as it would be the case for local disks dEpinat et al.l 
l2008alfbl) . This agreement is also true for some of the galaxies 
that we have classified as merging systems (see section 143] and 
individual comments in Appendix [A]). For mergers, the use of 
a fit relying on a disk hypothesis is undoubtly less meaningful 
than for relaxed systems. However, this fitting method enables 
to study the whole sample in an homogeneous way. 



lEpinat et alj (l2009h have shown that observation of high red- 
shift disk galaxies, strongly affected by beam smearing effects, 
have a disk inclination difficult to retrieve. The uncertainty on the 
inclination is the main source of errors in the rotational velocity 
determinations and thus on the Tully-Fisher relation, on the mass 
determination and on the mass assembly history. The extreme 
cases are provided by galaxies appearing face-on or edge-one 
from broad-band imagery: the radial velocity field of a face-on 
galaxy may show a velocity gradient indicating that the disk is 
not face-on; alternatively, an edge-on disk appears thicker and 
leads to an underestimation of the inclination. In order to check 
the disk inclinations determined by GALFIT, we have compared 
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Fig. 4. Comparison of Ha SFR with 02 (left) and photometric (right) SFRs. The dashed line indicates y - x. Black dots, red squares, 
green triangles correspond respectively to objects classified (see section 14.1b as rotating disks (RD), perturbed rotators (PR) and 
merging systems (MS). 




them to the ones computed using a simple correction provided in 
[Epinat et al. (2009). They modeled inclined thin disks having ex- 
ponential or flat truncated luminosity profiles with scale lengths 
from 2 to 6 kpc, observed at redshift z = 1 .7 with a seeing of 
0.5" and a sampling of 0.125". They introduced a parameter C 
in order to recover the true inclination / from the seeing FWHM 
(s) and the measured major (am) and minor (bm) axis lenghts 
(FWHMs) estimated in arcseconds from a 2D gaussian fitting : 



(7) 



where a and b are actual major and minor axis lengths (FWHMs) 
on the sky plane. The parameter C represents the fraction of 
the seeing by which both meas ured minor and m ajor axis are 
quadratically overestimated. In Epin at etaTI (|2009'), it is shown 
that C = 1 .014 ± 0.002 for an exponential luminosity profile and 
C - 1.015 + 0.010 for a flat and truncated luminosity profile. 
This means that these smooth luminosity profiles can be satisfac- 
torily described by a gaussian function. We have used C - 1.014 
to confirm the inclinations determined from the GALFIT fitting 
procedure. We have used the inclination determined from direct 
axis measurements for some objects for which GALFIT failed 
to converge into a realistic model. Except for those galaxies, 
GALFIT typically gives a disc component with a Sersic index 
of one and typical associated uncertainties ranging from 0.2 to 
0.5. 

In order to evaluate error bars on the inclination, the uncer- 
tainty on both major and minor axis lengths has been taken into 
account. Assuming that the error A is the same for both axes, 
we have computed upper and lower limits for the inclination re- 
ported in Table|4]following: 



b-A 



a + A 
b + A 



a - A 



(8) 



(9) 



The typical uncertainty is A = 0.25" (about one third the seeing). 

3.3. Kinematics of tiie ionized gas 
3.3.1 . Ha kinematics 

The SINFONI data were analysed using routines written under 
IDL environment 0. Our main goal was to determine the kine- 
matics of the ionised gas using the strongest emission line. Ha. 
In order to increase the signal to noise ratio without degrading 
the spatial resolution, we have performed a sub-seeing spatial 
gaussian smoothing (FWHM=2 pixels) on the data cube (the 
mean seeing is around 0.65" i.e. 5 pixels). For each galaxy of the 
sample we fit the spectrum around Ha to a single gaussian func- 
tion in order to characterize the emission line, and a pedestal to 
characterise the continuum. To minimize the effects of sky lines 
residuals on line parameters determination, the sky spectrum has 
been used in the fitting procedure in order to lower the weight of 
strong emitting sky lines wavelengths. From these fitting tech- 
niques it was possible, for each source, to recover the line flux 
map, the velocity field and the velocity dispersion map after cor- 
recting for instrumental dispersion (around 2.75A) measured on 
sky lines. 2D error maps have also been derived for each quan- 
tity from the fitting procedure. These are statistical errors that 
take into account the error spectrum and that indicate the accu- 
racy of the fit for each parameter. A signal-to-noise ratio map 
has been computed, the signal beeing the intensity of the fitted 
line and the noise beeing the standard deviation of the residual 
spectrum. These maps are presented in Figures [ATI to I A. 9l The 
velocity maps have been masked using several criteria: we have 
imposed (i) the line width to be larger than the one of the spectral 
PSF (the major part of the white pixels in the velocity dispersion 
error maps), (ii) the uncertainty on the velocity to be less than 
30 km s^^ and (iii) the signal-to-noise ratio to be larger than ~ 2 
(1.5 for the faintest object VVDS020182331). We also made a 
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Table 4. Kinematical parameters 



VVDS ID 


ZHa 


Inclination 


Position Angle 




V, 


Residuals {V) 


X 


V,/r, 












1 — 1 

km s 


km s 




km s kpc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


020116027 


1.53022 


A 4 4- 1 J 

44 

-23 


219 ± 8 


0.12 


32 


8 


2.4 


30 




l.ZZojZ 


52 _,3 


0/1/1 -L. A 


V.J J 


1 "3/1 


1 c 


1 




020147106 


1.51949 


37 

-37 

47 


307 ± 11 


0.12 


30 


3 


0.7 


29 


020261328 


1.52891 


192 ±5 


0.77 


195 


12 


2.2 


30 


220596913 


1.26619 


66 


253 ±5 


2.80 


325 


16 


3.5 


13 


220584167 


1.46588 


49!' 


178 ± 1 


1.24 


280 


15 


6.6 


26 


220544103 


1.39659 


61 i 


205 ±7 


6.73 


762 


17 


2.8 


13 


220015726 


1.29300 


184 ± 1 


0.12 


323 


16 


6.8 


309 


220014252 


1.31014 


48 


137 + 3 


0.12 


103 


13 


2.6 


99 



(1) Source VVDS identification number, (2) Geocentric SINFONI spectroscopic redshift deduced from model fitting, (3) Morphological 
inclination from CFHT /-band images, (4) Kinematical Position Angle of the major axis from North to East deduced from model fitting, (5) 
Characteristic radius of the kinematical model, (6) Characteristic velocity of the kinematical model, (7) Weighted standard deviation of the 
residual velocity map, (8) Reduced of the fit, (9) Inner velocity gradient derived from the model. 



careful visual inspection in the spectrum to remove outer regions 
validating the criteria but for which the line was associated with 
strong sky residuals (see Fi gures IA3] and IA!6] | . 

3.3.2. Model fitting 

To analyze the velocity fields of high redshift galaxies, four 
toy rotation curve models used in the l i terature (an expon en- 
tial disk by f orster Schreiber et alj|2006l: ICresci et al.ll2009L an 
isothermal disk b y Spano et all 2008 , an arc tangent function by 
Puech et al.l l2008l and a ilat rotation curve bv lWright et al.ll2007L 
2009 described hereafter) have been tested on a sub-sample of 
around 150 data cubes of nearby galaxies part of the GHASP 
samp le which have been projected at high redshift dEpinat et al.l 
|2009|) . These models consist of two parameters rotation curves 
that have various shapes, mainly the plateau (decreasing, in- 
creasing or flat) and the presence or not of an inner velocity 
bump. By comparing the p arameters de t ermin ed from high res- 
olution and projected data, Epinat et al.' ('2009') have shown that 
the simple model used by Wright et al. (2007, 2009) leads to the 
best estimation of the kinematical parameters in average. The 
high resolution rotation curve (velocity V as a function of the 
radius r) of this model is given by: 



rt 



when r < r, and 



(10) 



(11) 



when r > r,. 

In the present work, this fit was attempted for all sources on 
their velocity field using the associated velocity error map in or- 
der to evalutate which ones are compatible with rotating disks. 
For objects that we have classified as merging systems (section 
14. Il l, fitting rotating disks on substructures has been attempted 
and discussed in AppendixlAl The high resolution model is con- 
volved with the seeing measured on the PSF stars and by the 
2-pixels gaussian smoothing applied on the data cubes. The pa- 
rameters of the model are the position angle of the major axis 
(measured from North to East), the inclination, the center, the 
systemic velocity, the plateau of the rotation curve V, and the 
radius r, at which the plateau is reached. 



This model fi tting procedure described in details in 
lEpinat et alj (l2009l) is based on a minimization. One main 
concern for fitting the data is that the actual high resolution 
line flux distribution is unknown. We have assumed that the 
observed line flux is representative of the high resolution one. 
Furthermore, due to projection effects, kinematical models are 
affected by a degeneracy between the rotation velocity V{r) and 
the inclination. Even for local resolved galaxies, this dege neracy 
leads to large error bars on kinematical inclinations (Epi nat et alJ 
2008a, b). Due to the low sampling, this degeneracy can not be 
ruled out for high redshift galaxies. We have fixed the inclination 
to the morphological value determined from /-band continuum 
images as described in section l372l 

From the best fit model, we have computed a model veloc- 
ity dispersion map. This model velocity dispersion map only 
includes the enlargement of the profiles due to the unre- 
solved velocity gr adient because of beam smearing effects (see 
lEpinat et alj l2009l for analytical computation). It does not in- 
clude enlargement due to the physics of the gas. Thus by sub- 
tracting quadratically the model to the measured velocity disper- 
sion map, we are able to recover the local velocity dispersion 
noted hereafter ctq. In order to take into account the uncertainty 
on the velocity dispersion map while computing the mean veloc- 
ity dispersion, we have used a weight inversely proportional to 
the error. The same procedure has been applied while computing 
the mean velocity residual. 

The results of this fitting procedure are presented in Table 
|4] It contains the redshifts determined from SINFONI data, the 
morphological inclination used for the kinematical model, the 
kinematical position angle, the rotation curve parameters r, and 
Vi, the mean velocity residual, the reduced of the fit and the 
inner velocity gradient computed by the models. The center has 
been fixed to match the maximum of the flux distribution except 
for VVDS220584167 and VVDS220014252 (see Appendix|l| 
since the determination of the kinematical center is seve raly af- 
fected by the low spatial resolution as underlined in Epina t et aU 
(I2009I) . Since no clear turnover is reached in the velocity fields 
and due to a symmetrization process, the center could be off- 
set by more than 0.5"when let free, which is unrealistic. This 
biases the maximum velocity determination since the extent ap- 
pears to be larger on the one side than on the other. The bias 
depends on the offset. It also changes the value of the geocentric 
spectroscopic redshift to the value corresponding to the center. 
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VVDS ID 



(1) 



y 

' max 

km 
(2) 



134 

30 ^j; 

194 
177 
280 
146 
323 
103 +?? 



kpc 
(3) 



''1/2 

(4) 



km s^' 
(5) 



(6) 



-0.3 
+2.6 
-0.8 

-0.1 
+4.9 
-1.3 
+ 1.1 
-0.6 
+8.5 
-2.3 
+ 1.1 
-0.5 

-5.5 
+0.8 
-0.3 



10'" Mo 
(7) 



rrr 



-1.9 
+0.1 
-0.6 
+ 1.8 
-0.4 
+2.7 

+1-.1.9 

-5.0 

+2.3 

-3.1 

+5.7 

-1.6 

+ 1.0 

-3.6 



10'" Mo 
(8) 



4^1 

27.6!" 
19 6 

13-0 
21 1 

10.3 

16.7 !~ 3 
10.2 



Me 
10'" Mo 
(9) 



+0.6 
-0.1 
+ 1.4 
-0.6 

-0.1 
+9.7 
-1.7 
+ 1.4 
-0.7 

20.0 +^'32 

5 4 +'-'5 
•'■^ -0.7 

15 5 +°° 

2 3 +-■' 
^■•^ -0.6 



10'" Mo 
(10) 



(11) 



-0.0 
+8.4 
-1.1 

-0.0 
+ 1.9 



-Ml 
-11.6 
+1.4 
-0.5 



Mhaio 
10'" Mo 
(12) 



Class 



(13) 



020116027 
020182331 
020147106 
020261328 
220596913 
220584167 
220544103 
220015726 
220014252 



9.7 

6.3 

9.1 

6.4 

12.7 

10.9 

10.9 

6.4 

9.2 



2.9 
4.5 
2.0 
1.2 
8.7* 
7.9* 
5.5 
2.7 
4.4 



45 ± 13 
71 ± 34 
80 ±8 



55: 
76: 
47 : 
70: 
38: 
92: 



17 
20 
22 
17 
25 
19 



0.7 
1.9 
0.4 
3.5 
2.3 
5.9 
2.1 
8.5 
1.1 



1.2 
5.8 
1.7 
0.6 
8.5 
12.1 
5.1 
6.2 
6.1 



0.2 
2.6 
0.2 
5.6 
9.3 



5.3 
1.5 



+i:g 
-1.1 

27.4 



14.0 
3.7 
1.1 
4.9 
1.2 
7.9 



-5.7 
+ 10.1 
-7.4 
+2.3 
-1.7 
+1.3 
-0.8 
+2.7 
-2.1 
+2.1 
-1.0 
+3.7 
-3.0 



0.0 
1.8 
0.0 
0.4 
2.5 

18.1 
1.1 

13.2 



-11.6 
0.3 +"-^ 



0.5 
43.9 

0.4 
110.1 
98.7 
346.1 
51.0 
591.3 
19.4 



MS 
PR 
PR 
PR 
MS 
RD 
MS 
RD 
PR 



(1) Source VVDS identification number, (2) Maximum velocity deduced from model fitting, (3) Maximum radius of the kinematical maps, (4) 
Half light radius corrected for the beam smearing (an asterisk indicates that it has been computed from the /-band image), (5) "1/error"- weighted 
mean local velocity dispersion, (6) Ratio of the maximum rotation velocity (2) over the local velocity dispersion (5), (7) Stellar mass, (8) Total 
dynamical mass, (9) Rotation mass, (10) Dispersion mass, (11) Ratio of the dynamical masses, (12) Halo mass, (13) Kinematical classification 

(RD: Rotating Disk, PR: Perturbed Rotation, MS: Merging System). 



However, the position angle of the major axis remains a stable 
parameter. The centers used are displayed in Figures lATI to lA.gl 
The maximum velocity suggested by the models is sytematically 
higher than the maximum velocity deduced from the maps com- 
puted as 



2 X sin / 



(12) 



where v,„ax and v„ 



are the maximum and minimum values of 
the velocity field and / is the inclination (see Table |6]l. Indeed, 
the maximum velocity of a rotating disk is observed along the 
major axis. Due to beam smearing effects, the velocities mea- 
sured along the major axis are lowered by the contribution of 
off-axis regions. The beam smearing may also affect the veloc- 
ity gradients observed in merging systems velocity fields. For a 
given sampling, this effect is more pronounced for small galaxy. 

3.3.3. Dynamical masses 

To follow mass assembly accross the cosmic time, it is worth 
comparing dynamical masses and stellar masses even if uncer- 
tainties for both are rather large. 

Given our high resolution 3D data on rotating disk candi- 
dates, the most robust mass estimate is the enclosed mass es- 
timate. However a mass estimate requires several assumptions. 
First of all, the system has to be relaxed. The geometry of the 
system also has important consequences. If we consider that the 
mass is principally contained in spheroidal haloes, the enclosed 
mass Me at radius R can be expressed as follow: 



(13) 



V beeing the velocity at radius R. 

For some of the galaxies of our sample, there is no strong 
evidence for rotation and their kinematics are compatible with 
slowly rotating spherical systems, since the velocity gradient is 
lower than the velocity dispersion. For spherical, bound and dy- 
namically relaxed systems with random motions, the virial mass 
is more appropriate than the enclosed mass: 



virial 




log(Mstar) 



' 1 

Isunl 



12 



(14) 



Fig. 5. Comparison between total dynamical (y-axis) and stellar 
(x-axis) masses in log units. Black dots, red squares, green tri- 
angles correspond respectively to objects classified (see section 
14.11 ) as rotating disks (RD), perturbed rotators (PR) and merging 
systems (MS). The solid line represents y = x. 



The parameter C depends on the mass distribution and the ge- 
ometry of the system, ri/2 is the half light radius and cr is the 
mean random velocity. iBinnev & Tremain3 (l2008l) suggest that 
C = 2.25 is an average value of known galactic mass distribution 
models. 

Since our objects clearly show both rotation and dispersion 
motions, one might take into account both components to com- 
pute the total dynamical mass. To d o that, we have to a pply an 
"asymmetric drift correction" term (Meur er et alJI 19961) which 
involves radial gradients of the gas surface density, the gaseous 
velocity dispersion, and the disk scale height. To simplifiy, we 
assume that the disk scale height radial gradient is null as well 
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Fig. 6. Comparison between rotation (y-axis) and dispersion (x- 
axis) masses in log units. Black dots, red squares, green triangles 
correspond respectively to objects classified (see section [431 1 as 
rotating disks (RD), perturbed rotators (PR) and merging sys- 
tems (MS). The sohd hne represents y - x. 



as the dispersion one since we do not see clear radial depen- 
dency in our data. Thus, assuming that the gas is in dynamical 
equilibrium and that the gas velocity ellipsoid is isotropic, the 
total dynamical mass can be expressed as: 



Gh^ 



(15) 



where h is the gas surface density disk scale length when de- 
scribed by a gaussian function and cr is the local velocity disper- 
sion at radius R. We note the two terms of equation [15] respec- 
tively Mg and Mo- since they respectively refer to rotation and 
dispersion masses. 

The first term Mg is equal to the enclosed mass within radius 
R. The second term M,j is similar to the virial mass. Indeed, the 
gas surface density is well described by a gaussian function with 
ri/2 — h V2 In 2. The virial mass is thus equal to the asymmetric 
drift correction term at 



R^{CI2\n2f'^ ri/2 ~ 1.2 r^o 



(16) 



In Table |5] we have computed M^/,,,,, Mg and Ma- at Riast, 
radius of the last point for which a velocity measurement is 
done. Since the actual rotation velocity amplitude mainly de- 
pends on the true inclination of the system, we have used the 
uncertainties on the inclination to compute upper and lower lim- 
its for the rotation mass computed from equation [13] We have 
also used the maximum velocity Vmax derived from our kinemat- 
ical modeling within the radius Ruui (Yt when Riast > r,). Since 
we assumed the velocity dispersion to be constant, we used the 
"l/error"-weighted mean local velocity dispersion corrected for 
beam smearing effects ctq to compute the dispersion mass. The 
half light radius ri/2 has been determined for each galaxy from 
a two dimension gaussian fitting on line flux maps provided in 
Figures lA.ll to IA.9I The half light radius measured on /-band 



CFHT images is in average about 1.1 the Wa one. For two galax- 
ies (VVDS220584167 and VVDS220596913), since the fit on 
the Ha line flux maps did not converge, we thus used the half 
light radius measured on CFHT /-band images divided by 1.1. 
For each galaxy, ri/2 has been corrected for the beam smear- 
ing using a quadratical subtraction of half the seeing disk. All 
these parameters are stored in Table [5] together with the dynami- 
cal classification defined in section 14711 and two indicators of dy- 
namical support {Vmaxlo'a and Mg/Ma-) on which the classifica- 
tion partly relies. /?/„„ and ri/2 beeing approximatively linearly 
correlated, these two dynamical support indicators are directly 
linked since: 



^=21n2fc 

Ma- \ ri/2 



V,„ 



(17) 



The use of the mass estimates instead of Vmax/o'o ratio enables 
to know which is the main contribution. 

For galaxies classified as mergers (see section 14.1b . this 
method to compute the mass is more approximative than for ro- 
tating disks. Indeed, the rotation velocity and the extent mea- 
surements used to compute the rotation mass component Mg as 
well as the half light radius measurement used to compute the 
dispersion mass component M„- are derived using disk hypoth- 
esis. However, since it is difficult to extract unambiguously the 
components, this method gives a trend (see Appendix [Al for in- 
dividual details). 

As shown in Figure [5] the total dynamical mass is higher or 
equal than the stellar mass within the uncertainties. The differ- 
ence between dynamical and stellar masses could be even larger 
if for some galaxies the turnover in the rotation curve is not 
reached or if the extent of the stellar distribution is much larger 
than the gas one. This trend shows that in most of the cases, the 
gravitational potential is mainly due to gas and dark matter. 

Rotating disks (dynamical classes are defined in section [431 ) 
have both large dynamical and stellar masses (around 10" M©) 
while mergers and perturbed rotators present a larger mass 
range. The comparison between dynamical and stellar masses 
is discussed for each galaxy in Appendix[Al 

In Figure [6] comparing rotation and dispersion masses, 
rotation-dominated systems are clearly distinguished from 
dispersion-dominated ones. Indeed, on the one hand, for the 
two galaxies classified as rotating disks (VVDS220584167 and 
VVDS2200 15726), the rotation mass is higher than the dis- 
persion mass by a factor larger than 13 showing that their 
gravitational support is rotation dominated. On the other hand, 
VVDS0201 16027 and VVDS020147106 are clearly dominated 
by random motions (Ma/Mg > 10) even if inclination effects 
may lead to slightly higher values. The five remaining galaxies 
with Mg/Ma ~ 1 show both rotation and dispersion gravitational 
supports. 

3.3.4. Halo masses 

We have estimated the total halo mass (M/m/q) using a spher- 
ical virialized collapse model (Peebles 1980; IWhite & Fr enkI 
ll991l:lMo & Whitell2002l) . We assumed that the plateau has been 
reached in our observations and that the plateau velocity Vmax 
traces the modeled halo circular velocity. 



Mhaio = 0.1 G ' (l+z)-'-' Vi 



(18) 



where Hq is the local Hubble constant, Q.„, is the local matter 
density, G is the universal gravitational constant and z is the red- 
shift of the galaxy. 
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Fig. 7. VIMOS optical spectra from which the original redshifts were derived (left) and SINFONI (right) ID spectra around the 
Ha line for the four VVDS-02h field galaxies. SINFONI spectra have been produced summing spectra of all spatial resolution 
elements with detected Ha emission. The main emission and absorption lines are marked and labeled for each sources. The line 
wavelengths are computed using the redshifts deduced from SINFONI data (Table|4]i. 



In the cases where our models indicate that the plateau is not 
reached within the observation, we have used the maximum ve- 
locity derived from our kinematical modeling within the radius 
Riasi instead of the plateau velocity of the model V, since the 
plateau is not constrained correctly. The values are available in 
Table 121 



If we assume that the length of the halo is the size that a 
galaxy should have in order that its enclosed mass is equal to the 



halo mass, then for our sample, the halo length is 113 
larger for fast rotating disks since: 



R 



halo 



= 0.1//o'^«''(l+^) 



-1.5 



: 73. It is 



(19) 



These total halo masses have no physical meaning for galax- 
ies which are not gravitationally supported by rotation. For 
galaxies having a ratio Vmax/cro larger than unity, the total 
halo masses range from 0.2 to 6.0x10'^ M©. These values are 
much larger than the measured dynamical masses. At redshift 
z = 1.5, a galaxy with the same maximum rotation velocity 
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Fig, 8. Similar as Figure|7]for the five VVDS-22h field galaxies. 



than a z = galaxy is supposed to have a total halo mass 
lower by a factor ~ 4. Our rotation dominated sources have 
Mo/Mi,aio ~ 0.049 + 0.015 in average. This value cannot be di- 
rectly compared to the fraction of baryonic matter to total matter 
Q.hlQ.,„ - 0.171 (K omatsu et a l. 2009) because Mg may con- 
tain a noticeable fraction of dark matter. Indeed, for comparison, 
typical halo masses measured from mass model derived from 
rotation curves at the last HI radius for nearby galaxies larger 
than ISkpc is around 25 x 10'° Mq ( Blais-Ouellette et al.ll200ll: 
ISpano et al.l2008.) . Taking into account that the values computed 
for local galaxies are measured, in average, at a radius twice as 
large as the sample of high redshift galaxies, the total cosmo- 
logical halo mass for high redshift galaxies are still much larger 



than the total halo mass measured for local galaxies from rota- 
tion curves. 

3.4. Search for type 2 AGN and metallicity 

Along with the VIMOS optical spectra, we show in Figures [T] 
and |8] the integrated ID SINFONI spectra around the Ha line, 
resulting from the sum of all spatial resolution elements in the 
datacubes at the position of the sources. VIMOS and SINFONI 
combined spectra of the nine observed galaxies are shown in 
Figure|9] These are median spectra that have been computed us- 
ing respectively [On] flux and Ha flux to normalize each in- 
dividual spectrum. One can clearly identify [Nii] /16548 and 
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Fig. 9. Combined spectrum of the nine VVDS galaxies at 1 .2 < 
z < 1.6 observed with VIMOS (top) and SINFONI (bottom), in 
arbitrary flux units and rest frame wavelengths. The position of 
the main emission and absorption lines are indicated in red and 
blue respectively. 



[Nil] /16584 emission lines, as well as the [S ii] doublet at 6717 
and 6731 A in the SINFONI combined spectrum. The low [Nii] 
A6584/Ha and [S ii] /1/167 17,673 l/Hff ratios (=0.13 and 0.23 re- 
spectively) indicate the lack of AGN contribution in our sample. 
We further performed a careful inspection of the spatial distri- 
bution of the [N ii] ^6584/Hq' ratio for each galaxy. [N ii] /16584 
line flux maps were derived similarly to Ha ones (see section 
13.3.1b . However, since the [Nii] /16584 flux is rather weak, we 
used the Ha line to constrain the line position. We do not find any 
significant peak in the flux ratio maps, contrary to Wrig ht et al.l 
(2009) who found an AGN in two of their six z ~ 1.6 galax- 
ies. This lack of AGN contribution is further supported by other 
indicators. 

- the quite low mean value of reddening (E{B - V) ~ 0.3 or 
Ay ~ I) measured in our sample galaxies is against any dusty 
AGN alternative 

- we do not detect any strong [Neiii] /13869 or [Nev] /i3426 
emission in the VIMOS composite spectrum which would be 
a clea r signature of Sey2-like activity (e.g. IZakamska et al.l 
I2003h 

- no evidence for Sey2-like activity is found in X-ray data. 
Indeed, we have checked in the HEASARC archive that no 
X-ray counterpart was yet detected at less than 30" from our 
objects. Moreover, when existing, the counterparts are closer 
to brighter objects. 

Note also that type 1 AGN are excluded from the parent VVDS 
sample on the basis of the presence of broad emission lines in 
the VIMOS spectrum. 

The [N ii] /16584/Ha ratio can be used to estimate roughly the 
typical metallicity of our sample. Using the last calibrations pro- 
posed bv lPerez-Montero & ContiM (l2009l) . we have estimated a 




Fig. 10. Maximum rotation velocity V,„a.x deduced from rotating 
disk model fitting as a function of "l/error"-weighted mean lo- 
cal velocity dispersion corrected for beam smearing effects ctq. 
Black dots, red squares, green triangles correspond respectively 
to objects classified as rotating disks (RD), perturbed rotators 
(PR) and merging systems (MS). The soUd line represents y - x. 

rather low oxygen abundance l2+l og(0/H) - 8.37, cor respond- 
ing to half the solar metallicity (see lOuevrel et al.l2009l for a full 
discussion on the mass-metallicity relation). 

4. Results 

4.1. Kinematical classification 

It has already been observed by other authors studying 

-spectroscopy jLaw et al.1 



z > 1 star-forming galaxi e s with 3D-spec 
20071 IW right et al.' ^2007'. l2009t iGenzel 



et al.l 120061 120081: 



Forster S chreiber et al. 2006 yCresci et al.l2009h that the gaseous 



velocity dispersion of these objects is larger than at lower red- 
shifts (Yang et al., 2008). As shown from kinematical model- 
ing, an unresolved velocity gradient cannot a ccount for high 
velocity dispersion values (lEpinat et al.l 20091). Th e kinemat- 
ical classification as performed by iFlores et"an (|2006|) on 
FLAMES/GIRAFFE data relies mainly on the ability to observe 
a signature of unresolved velocity gradient in the velocity disper- 
sion map: rotating disks are defined as systems showing a regu- 
lar velocity gradient and a central peak in the velocity dispersion 
map, where the actual velocity gradient is the steepest; systems 
with an off'set velocity dispersion peak from the center, or with- 
out peak, are defined as perturbed rotators; the galaxies with no 
regular velocity gradient and no velocity dispersion central peak 
are classified as galaxies having complex kinematics. For sys- 
tems with a high local velocity dispersion of the ionized gas, the 
velocity dispersion argument is thus not valid anymore. In Figure 
nn we show examples of rotating disk models for two objects 
from our sample illustrating this point. For VVDS2200 15726 
(first fine) and VVDS2200 14252 (second line), the model ve- 
locity fields (first column) show that the kinematics of these ob- 
jects reasonably match rotating disks (the second column shows 
the residual velocity fields). For VVDS220015726, the signature 
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Fig. 11. From left to right: Model velocity field, velocity field residuals, model velocity dispersion map and velocity dispersion 
map quadratic residuals. North is up. East is left. The center of the models is marked on each kinematical map with a black and 
white double cross. The solid line represents the orientation of the kinematical major axis and ends at half light radius measured on 
I-band images. The flux contours are overplotted. 



in the velocity dispersion map (third column) correctly matches 
the unresolved velocity gradient signature since the velocity dis- 
persion map central signature (see Figure [ATZt is well removed 
when subtracting quadratic ally the model (last column of Figure 
fTTT i. For VVDS220014252 which has a higher velocity disper- 
sion, the model velocity dispersion map absolutely does not ac- 
count for the observed velocity dispersion as demonstrated by 
the resemblance between corrected (last column of Figure fTTT i 
and uncorrected (see Figure |AT2] i velocity dispersion maps. 

The ~0.65" seeing continuum images from CFHT do not 
all ow a morphol ogical classification as accurate as performed 
by iNeichel etaTI ((2008) from HST data on IMAGES sample. 
However, relying on line flux morphology, the kinematical maps 
derived from our SINFONI data and from attempts to fit a ro- 
tating disk model, we are able to derive a first order kinematical 
classification. 

We defined rotating disks (RDs) as systems for which the 
rotation mass is higher than the dispersion mass (rotation dom- 
inated systems, see the two parameters Ymaxicro and MelM^ 
linked by equation [TTl and given in Table [Sj as shown in Figure 
|6] and for which both velocity field and velocity dispersion map 
are described by rotating disk models. Mergers systems (MSs) 
are defined as systems with two spatially resolved components 
in the Ha map and for which the velocity field is not described 
accurately by rotating disk models. Other systems are classified 
as perturbed rotators (PRs). They mainly show high velocity dis- 
persions, higher than around 60 km i and their velocity fields 
show some peculiarities not described by rotating disks. 

Figure [TO] shows a correlation between the maximum rota- 
tion velocity y,„„v and the "l/error"-weighted mean local veloc- 
ity dispersion corrected for beam smearing effects ctq, in partic- 
ular for rotating disks and perturbed rotators. The most massive 
galaxies (see Table |5j are our two RDs. This would imply that 
the most massive disks are stable earlier or that stable disks can 
already be formed at these redshifts. Random motions convert in 
organized motions (rotation) more rapidly for massive systems. 



Table 6. Magnitudes and maximum velocities 



VVDS ID 










^iKis.v, model 


^max, map 




mag 




mag 


km 


km 


(1) 


(2) 




(3) 




(4) 


(5) 


020116027 


-22.7 ± 


0.4 


-21.8 ± 


0.2 


32 +f 


23 !f 


020182331 


-22.3 ± 


0.2 


-21.7 ± 


0.1 


134 




020147106 


-22.8 ± 


0.1 


-22.1 ± 


0.1 


30 


21 « 


020261328 


-21.6 ± 


1.0 


-21.1 ± 


0.3 


194 


82 


220596913 


-23.2 ± 


0.5 


-22.6 ± 


0.1 


177 


153!' 


220584167 
220544103 


-24.2 ± 
-23.2 ± 


0.5 
0.5 


-23.1 ± 
-22.5 ± 


0.1 
0.2 


280 +°1 
146 


198 % 
108 *f 


220015726 


-23.2 ± 


0.5 


-22.1 ± 


0.1 


323 


189!" 


220014252 


-23.3 ± 


0.5 


-22.6 ± 


0.2 


103 


84 !5? 



(1) Source VVDS identification number, (2) AB absolute magnitude in 
rest frame AT-band, (3) AB absolute magnitude in rest frame i?-band, 
(4) Maximum velocity derived from the model, (5) Maximum velocity 
derived from the maps (equationll2t. 



From our classification, we obtain three clear major merging 
systems, two clear rotating disks and four perturbed rotators that 
could be either perturbed by minor merging events or by contin- 
uous gas accretion. The classification is summarized in column 
(10) of Tableland individual detailed comments on the kine- 
matics and dynamics of each galaxy are provided in Appendix 



El 



4.2. The Tully-Fisher relation 

The Tully-Fisher relation (hereafter TFR) jTuUv & FisheJI 19771) 
links the maximum rotation velocity to the actual flux or mass 
of a galaxy (via their logarithm). Extensive work has been pro- 
duced in order to probe the evolution of the TFR with cosmic 
time in both Z?-band and /T-band up to redshifts z ~ 1.4 (e.g. 



B. Epinat et al.: Integral field spectroscopy with SINFONI of VVDS galaxies 



15 




1 00 

[km s-'] 



1 000 




1 00 
V„.o. [km s-'] 



1 000 




1 00 
V„.ax [km s-^] 



1000 



Fig. 12. A'-band (top), B-band (middle) and stellar mass (bot- 
tom) TFRs for our nine z ~ 1 .4 galaxies. Black dots, red squares, 
green triangles correspond respectively to objects classified (see 
section l4m i as rotating disks (RD), perturbed rotators (PR) and 
merging systems (MS). The maximum velocity has been deter- 
mined from the rotating disk fit modeling. Top: th e dotted line 
repre sents the local reference K-hand TFR used bv iPuech et al.l 
(l2008t) derived from Ha on a SDSS sub-sample. Top and mid- 
dle: the solid line represents respectively the local /T-band and 
B-ban d TFRs derived from a local HI sample by iMever et al.l 
(l2008h whereas the dashed lines represent the dispersion of the 
HI galaxies around these local relati ons. Bottom: the sol id line 
represents the z = TFR derived bv iBell & de Jond (12001 ). The 
da shed Unes represent the z = 2.2 TFR derived from SINS data 
bv lCresci et all (l2009l) . 



[Fernandez Lorenzo et al.ll2009l for references). Since it requires 
a large sample of galaxies (more than ~50 galaxies), high red- 
shift TFR were derived fro m long-slit sp e ctrosc opy samples, ex- 
cept th e ones derived by (i) iFlores et al.l (l2006l) and lPuech et alj 
(120081) from FLAM ES/GIRAFFE IFUs data at z ~ 0.6; (ii) 
Cresci et al.1 (|2009|) fr om SINFONI data at z ~ 2.2 and (iii) 
Swin bank et al.l ( 20061) from lensed galaxies at z ~ 1 observed 



with GMOS. 

In the frame of the present work, due to the small number 
of objects, we only look at the distribution of the different kine- 
matical types accross the local TFR. We pre sent the local TFR 
derived from HI data by Mever et al.! (l2008b (solid line) since 
they provide both jiT- band (Figure [12] top ) and B-band (Figure 
[T2I middle) rel ations. iPuech et al. (l2008l) used a subsample of 
SDSS sample ('Pizagn o et al.l 1200 7) to derive a local reference 
TFR. The latter has been derived using Vgo velocities instead 
of maximum velocities (from Ha kinematics) since they argue 
that their maximu m v elocity is statistical ly not measured (see 
iPuechet all [20081 and lPizagno etal.ll2007l for details). We thus 
also plot this /T-band reference TFR in Figure [12] top (dotted 
line). 

For local samples with Ha kinematical data, usually one has 
to exclude galaxies with low inclination because the uncertainty 
on their maximum velocity is high, but also galaxies with high 
inclination because of dust extinction in the disk that can bias 
the maximum velocity determination. For high redshift, as men- 
tioned in section 13.21 the inclination is poorly known, thus we 
expect to observe both face-on and edge-on systems that we are 
not able to identify. 

Absolute magnitudes in the desired rest-frame band (using 
the appropriate filter response curve) and stellar masses with 
their associated error bars for our objects were obtained from 
the same spectral energy distribution modeling than described 
in section[TT]using the GOSSIP software (Franzetti et al. 20081). 
Stellar masses and the associated errors are provided in Table 
[5] whereas rest-frame B- and /T-band absolute magnitudes and 
their associated errors determined from the PDF are provided 
in Table [6] Uncertainties on ^T-band magnitudes are larger than 
for B-band. They result from the lack of any measurement of K 
rest-frame (~ 3.5//m) to constrain efficiently the SED fitting. 

In Figures [12] top and middle, we respectively show the K- 
and B-band TFR for our nine z ~ 1.4 objects. The maximum 
velocities V„,ax have been derived from the rotating disk mod- 
els. Vmax is systematically higher when coming from the models 
that take into account the beam smearing eff'ects (see Table |6]l. 
On these Figures, to distinguish the three kinematical types we 
affected different colors and symbols (see caption). We see that 
the RD galaxies follow accurately the local B-band TFR as well 
as the PR with the highest y„„„. These galaxies show a better 
agreement with the SDSS /T-band TFR rather than the HI one. 
The two objects with no consequent rotation observed are the 
furthest from the local TFR (both K- and B-band). For the one 
classified as PR, this can be accounted for by the inclination. As 
for local galaxies, these galaxies have to be excluded from the 
TFR analysis because of the uncertainty on the inclination. The 
other objects classified as PR or MS stand on a restricted area on 
the TFR and present some correlation (if we exclude the slowest 
rotator) parallel to the local TFR, both in B- and /iT-bands. The 
three ki nematical types see m to behave as the three kinematical 
types of lPuechetaLl (120081) . 

Given the strong evolution in luminosity with redshift and 
the fact that the B-band rest-frame magnitude is sensitive to short 
bursts of star formation, we also consider in Figure [12] (bottom) 
the stellar mass TFR which has the advantage of probing more 
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Fig. 13. Ratio of the maximum rotation velocity Vmax (from 
modeling) with the mean local velocity dispersion corrected for 
beam smearing effects ctq as a function of the maximum rotation 
velocity. Black dots are our rotating disks (RD) and perturbed 
rotators (PR). The grey cone corresponds to the area cove red by 
local rotating disk galaxies (GHASP sample. lEpinat et ar]|2009l) . 
Galaxies fro m the literature are also plotted: dark green trian- 
gles are fro m Law et all (l2007l 120091) . light green star is from 
|Sto-ketal.'(2008'), pink up side down triangles are rotating disks 
from Wright et al. (2007, 20 Q9) and red s quare s are rotating 
disks from the SINS sample iGenzel et alj dioOSl) : ICresci et al] 
(1200 9). The maximum rotation velociti es and the mean local ve- 
locity dispersion for the galaxies from iLaw et al. I (I2007ll2009t) 
have bee n corrected f or beam smearing using the recipes pre- 
sented in lEpinat et alj (|2009), and a mean statistic inclination of 
45° has been used to deproject the velocities in the galaxy plane. 



closely the actual stellar mass build-up of galaxies. As for the 
absolute magnitude T FR, the two rotating disks follow the stellar 
mass TFR traced by ICresci et a n (120091) between z ~ 1.5 and 
z ~ 2.5. 

A larger sample will be needed in order to provide more ad- 
vanced conclusions on TFR evolution. 



5. Discussion and conclusions 

Using high spatial resolution near infrared integral field spec- 
troscopy, we have produced a detailed analysis of the kinemat- 
ical maps of nine galaxies with 1.2 < z < 1.6 selected from 
the VVDS. We have identified three major mergers systems, two 
massive rotating disks with total dynamical mass > 10" M©, 
and four galaxies with an indication of rotation but with per- 
turbed velocity maps, two of them with h igh gaseous ve locity 
dispersion unseen in local isolated disks (Epinat et al."2009). Our 
sample, although still small but comparable in size to other sam- 
ples in the literature, has the advantage to be selected from an un- 
biased sample drawn from the /-band magnitude selected VVDS 
and with [Oil] 3727 EW > 50 A. 

The effects of beam smearing on the data (mean seeing 
around 0.65" corresponding to around 5 kpc at these redshifts) 



have been taken into account in the interpretation of the kinemat- 
ical maps in particular to correct the velocity dispersion maps 
from unresolve d velocity gradients, us i ng mo del fitting meth- 
ods as done by iForster Schreiber et alj (l2006l) or IWright et~al] 
(I2007h . 

Finding three major merging systems, representing 30% of 
our sample, is a high occurence. We can co mpare this pair 
fractio n to the measurements at lower redshifts: Ide Ravel et al.l 
(l2009h find that at z ~ 1 about 10% of galaxies with a stel- 
lar mass M» > 10"^ Mq are in pairs likely to lead to a major 
merger, making the merger rate in our higher redshift sample 
significantly higher When taking into account typical merger 
timescales of ~ 0.5 - 1 Gyr, the high rate in our sample means 
that over a ~ 3 Gyr period virtually all galaxies with stellar 
mass > 10'" Mq will have undergone at least one major merger 
event. Although the observed sample is small, this result is com- 
parable to the frac tion of mergers in o ther samples at sim ilar 
(IWright et al]|2009i) or higher redshifts dShapiro et alj|200 8). In 
two of our merger systems, we can identify that one of the galax- 
ies is compatible with a rotating disk, indicating that disk galax- 
ies are indeed involved in merging events at these high redshifts. 

The two massive rotating disks, with velocity dispersion 
maps indicating stable disks, confirm that massive rotating disks 
with total dynamical mass M > 10" Mq are already well assem- 
bled at z ~ 1.5. The stability is underlined by the agreement of 
these galaxies with local Tully-Fisher relations obtained from lo- 
cal samples. These disks obviously need to have been assembled 
at earlier cosmic times, therefore requiring an early process of 
formation. These objects are compara ble to some o f the m assive 
disks observed at s imilar redshifts by Wright et al. ( 2009 ), or at 
higher red shifts by iForster Schreiber et al. (2006). Genzel et al.l 
(l2008h and lLaw et al. (2009). However, we note that these disks 
are only about 20% of our sample. 

Perturbed rotators are the dominant population in our sam- 
ple with four galaxies over nine. These objects show perturbed 
velocity maps, two of them with high velocity dispersion, with a 
y/cry lower than for typical rotating disks (Figure[T3j, indicating 
that the gravitational support is not mainly due to rotation. These 
objects are clearly a separate class not observed at lower red- 
shifts. Their dyna mical state is comparabl e to som e of the galax- 
ies ob ser ved bv IForster Schreiber etal] (|2006|) . IWright et~a!] 
(120091) . or lLaw et al.1 (l2009h . At a resolution of a few kilopar- 
secs, we argue that it is difficult to disentangle whether the high 
velocity dispersion is produced by instabilities resulting from 
gas accretion e.g. along the filaments of the cosmic web, rep- 
resenting early-stage disks with a high cold gas fraction which 
fragmented under self-gravity and collapsed to form a starburst 
llmmeli et al. 2004b), or by gas rich minor merger events (e.g. 
[Semelin & Combe s 2002). Indeed, both scenarios of continuous 
accretion of cold gas coming from cosmological filaments and 
frequent mergers of minor cold gas satellites may fuel the disk 
in fresh gas. If minor mergers (ranging from 10:1 to 50:1) are 
relatively significant, the relaxation processes will expulse the 
pre-existing stars from the disk to a spheroid structure or to a 
very thick disk not directly observable for high redshift galaxies. 
On the other hand, the formation of these structures would not 
perturb dramatically the kinematics of the disk, and thus would 
not be easily observable. Indirect signatures of the existence of 
a spheroid could be the stabilization of the gaseous disk, reduc- 
ing the star formation out of merger phases and thus reducing 
the amount of sub-structures like Ha or UV clumps in the disk. 
On the contrary, a very smooth accretion of diffuse gas coming 
from an extended halo or from the cosmological filaments will 
not be so efficient to form a stellar spheroid. Thus instabilities 
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would be more numerous in the disks and possibly observable 
through deep imaging addressing the formation of clumps. For 
very minors mergers (e.g. mass ratio > 100:1), the dwarf galax- 
ies will be dislocated by the tidal field once they experience the 
gravitational field of the main galaxy. Thus, torn by tidal field, 
the accretion of very small companions will resemble very much 
to diffuse gas accretion. For a baryonic mass galaxy of 10'" M©, 
these satellites correspond to mass lower than 10*^ Mq. If they 
exist, due to the lack of spatial resolution, these galaxies are not 
observed, neither in observations nor in numerical simulations. 

For two of the perturbed rotators observed in our sample 
(VVDS020147106 and VVDS220014252), the velocity maps 
are asymmetric and consistent with a small companion in the 
process of merging. More generally, these perturbed rotators in- 
dicate that disk build-up is fully on-going at the peak of cosmic 
star formation history. 

Star formation rates in our sample are high, with an average 
of ~ 55 Mq yr"', meaning that over the ~1 Gyr between z = 1.6 
and 1.2 more than 5 x 10'" M© would be formed if this rate 
is sustained, comparable to the average dynamical mass. This 
seems to indicate that any external accretion should last only a 
limited time. 

We have not found any evidence for type 2 AGN in our sam- 
ple, (type 1 were already excluded from our parent galaxy sam- 
ple). The average [Nii]/Ha' and [Sii]/Ha' ratios are typical of 
star-forming galaxies (Stasiriskaet al. 2006), and the [Nii]/Ha' 
2D maps do not show a central excess of [N ii] which could in- 
dicate an AGN. Although we cannot exclude low luminosity or 
highly obscured AGN, this seems to indicate that the AGN and 
star formation phenomena may be only weakly connected in the 
main star-forming galaxy population. 

The question of the assembly of galaxies via major dissi- 
pative mergers or internal secular processes has been recently 
highly debated in the literature. Genzel et al. (2008) advocated a 
secular process of assembly t o fo rm bulges and disks in massive 
galaxies at z ~ 2. iLaw et alj (12009.) concluded that the high ve- 
locity dispersions they observe in most galaxies at z ~ 2 may be 
neither a 'merger' nor a 'disk', but rather the result of instabil- 
ities related to cold gas accretion becoming dynamically domi- 
nant. Our data seem rather to indicate that several processes are 
acting at these epochs. Among them, merging seems to play a 
key role. On the one hand, we find a high 30% rate of close pairs 
of galaxies expected to merge in less than 1 Gyr, indicating that 
the hierarchical build up of galaxies at the peak of star forma- 
tion is fully in progress. On the other hand, in our sample the 
dominant 'perturbed rotators' may include a significant fraction 
of galaxies with minor mergers in progress or cold gas accre- 
tion along streams of the cosmic web, producing a high velocity 
dispersion. Minor merging may be an important phenomenon in 
the build-up of galaxies, which will require higher spatial reso- 
lution than currently available on 8-lOm telescopes, even with 
adaptive optics, to be confirmed. In large numerical simulations, 
major and minor mergers are indeed playing an important role 
( Bournaud et al.ll2b05i). in parallel to continu ous gas accretion 
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Semehn & Combesll2005l: iDekel et al.ll2009l) . Whether the two 



massive disks that we observe in our sample are the result of sec- 
ular processes with continuous accretion (Genzel et al. 2008) or 
of anterior hierarchical build up by major and/or minor mergers 
cannot be assessed as both processes take a relatively short time 
(less than a 1 Gyr) to complete. These arguments will need to be 
revisited from larger representative samples. 

In conclusion, the statistics of this representative sample, 
based on nine galaxies, show that there does not seem to be 
one single process driving the mass assembly in galaxies. Major 



mergers certainly play an important role, while the contribution 
of minor mergers is likely but will remain difficult to confirm. 
In addition, secular evolution with accretion which drives gas 
and stars in the central regions of galaxies remains a possibility 
to assemble bulges and disks early in the life of the universe. 
Furthermore, the absence of AGN type 2 in our sample indicates 
that the AGN phenomenon in high redshift galaxies is at best a 
short lived event. With the small samples observed to date, we 
cannot exclude that in the life of galaxies these processes will 
act at different epochs. 

Whether the differences between our sample and other exist- 
ing samples are mainly due to different selection functions will 
need to be investigated with larger well controlled samples. With 
the on-going MASSIV program at the VLT, we will be able to 
investigate these issues in more details. 
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Appendix A: Properties of individual galaxies 

In this appendix, we present the morphological, kinematical 
and dynamical properties for each galaxy derived from /-band 
images. Ha flux maps, radial velocity maps, velocity dispersion 
maps and their associated error maps (Figures I A. 11 to |A.9t . 
On each velocity map, the kinematical center used for the 
kinematical analysis is plotted with a black and white double 
cross. The kinematical position angle is also displayed and 
ends at half light radius ri/2 (Table |5]). The flux contours are 
overplotted on both velocity fields and velocity dispersion maps. 
In all figures. North is up and East is to the left. The kinematical 
parameters and the masses are stored in Tables|4]and|5] On error 
maps, pixels with values out of the scale range are displayed 
using white. 

VVDS0201 16027 

This galaxy shows a faint and irregular velocity gradient. 
The Ha flux map, the velocity field and the velocity dispersion 
map suggest two components separated of about 6 kpc (0.75")- 
Indeed, the velocities decrease from South to North accross 
the main component but decrease of around 15 km i"' when 
reaching the faint component. The separation between the two 
components is characterized by a higher velocity disp ersion 
that could indicate a merging event dAmram et alj 120071) or an 
unresolved gap in velocity between the two components. The 
velocity separation of the centers of these two components is 
close to km s^^. The me rging timescale would be around 
0.4 Gyr (from lKitzbichler & White 2008) . The faint component 
is not seen in the /-band CFHTLS image, but the system 
appears as distorted and elongated, leading to an inclination 
of the whole system of 50°. The very low maximum rotation 
velocity (50 km s~^) may indicate that the object is almost not 
rotating {Vma.x/crQ = 1.1) or that the main component is seen 
nearly face-on. Fitting a rotating disk model for this system 
has been attempted both with or without masking the northern 
component. In both cases, a low is obtained but is due to 
the low velocities measured in this galaxy. The cr-peak around 
70 km i"' located at 0.3" at the North-East of the maximum of 
the Ha flux is not accounted for by these models. Due to both 
the low velocity dispersion and maximum rotation velocity, this 
object is one of the less massive of our sample (5.5 x 10'" M©). 
This object is classified as a major merger, presumably in a 
premerging state since the SFR//„ is rather faint, indicating that 
the starburst event would not be yet engaged. 

VVDS020182331 

The flux of this galaxy is mostly concentrated in the central 
parts. This galaxy shows a perturbed velocity field. A strong 
sky line residual (14605 A) perturbs the line measurements and 
is responsible for both northern velocity dispersion and line 
flux peaks observed towards the North. The uncertainty on the 
velocity measurements is lower than 20 km i"' for the central 
arcsecond only (in diameter). In this central part, a rotating 
disk model is acceptable. The corrected velocity dispersion 
of this object is high (71 km s^^). We do not detect in Ha 
the close-by object revealed in CFHTLS best seeing images, 
indicating that this object is a foreground or background galaxy 
or that it is a companion with no star formation with which 
VVDS020182331 may interact. The latter explanation could 
explain the perturbed rotation observed in VVDS020182331. 
This object is slightly more rotation supported than dispersion 
supported (Vmax/o'o - 1.9) and has a dynamical mass of 
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Fig.A.l. VVDS0201 16027 maps. From left to right: /-band CFHT legacy survey best seeing image. Ha flux (top) and signal- 
to-noise (bottom) maps. Ha velocity field (top) and associated errors (bottom). Ha velocity dispersion map and associated errors 
(bottom) obtained from Gaussian fits to the SINFONI data cubes after smoothing spatially with a two-dimensional Gaussian of 
FWHM - 2 pixels. Except the instrumental spectral PSF, no correction has been applied to compute the velocity dispersion map. 
Kinematical maps have been masked using the following criteria: (i) the line width must be larger than the one of the spectral PSF 
(the majority white pixels in the velocity dispersion error map), (ii) the uncertainty on the velocity must be less than 30 km 5"' 
and (iii) the signal-to-noise ratio larger than ~ 2. The /-band image and Ha maps are color-coded with a linear scaling such that 
the values increase from light to dark. The integrated Ha flux is quoted in each Ha flux map in 10 '^ergs s ' cm"^. In error maps, 
values larger than the scale maximum are displayed using white. The angular size is indicated and the final seeing FWHM (including 
observational seeing and gaussian smoothing) is plotted as a circle on the velocity field. North is up. East is left. The center used 
for kinematical analysis is marked on each kinematical map with a black and white double cross. The solid line represents the 
orientation of the kinematical major axis and ends at half light radius measured on I-band images. The flux contours are overplotted 
on the velocity field, the velocity dispersion map and their associated error maps. The redshift derived from SINFONI data and the 
dynamical classification are quoted (RD for rotating disks, PR for perturbed rotators, MS for merging systems). 



1" 




Fig. A.2. Similar as Figure lATI for VVDS020182331. The signal-to-noise threshold for this galaxy has been set to 1.5. 



4.1 X 10 Mq. It is classified as a perturbed rotator. 
VVDS020147106 

The Ha flux is concentrated in the center. The size of the 
galaxy is twice the seeing (see Table [Til. This galaxy shows a 
faint velocity shear compatible with a face-on rotating disk. 



The velocity dispersion map does not show any peak but is 
higher on the eastern side (of about 20 km s '). The northern 
edge shows higher velocities of about 30 km s However the 
one pixel crown of the edges has higher uncertainties due to a 
low signal to noise ratio (lower than 4) and was masked to fit a 
rotating disk model. The fit indicates that a low velocity plateau 
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Fig. A.3. Similar as Figure|Aj]for VVDS020147106. 
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Fig. A.4. Similar as Figure |AT|for VVDS020261328. The East feature in the signal-to-noise map correspond to a badly subtracted 
cosmic ray. 



(30 km i"') is reached close to the center This low maximum 
rotation velocity suggests that the inclination of this system is 
even lower than suggested from the /-band image or that this 
object is strongly dispersion dominated {Vmax/cro - 0.4). Due 
to compactness of this object (ri /2 = 1.2 kpc) and compared to 
the total extent of Ha emission {Riasi = 6.4 kpc), the dispersion 
mass is estimated to 27.4 x lO'^M©, which is 15 times larger 
than the stellar mass (1.7 x 10'" M©). This suggests that this 
star-forming galaxy is embedded in a large dark matter halo 
or that it contains a large amount of gas. Even if it can not be 
excluded that this galaxy is a dispersion dominated spheroid , 
this galaxy is classified as a perturbed rotating disk since it has 
a high velocity dispersion (80 km .?"') and shows rotation in 
its perturbed velocity field. The asymmetries observed in the 
velocity maps could be accounted for by a minor merging event. 



VVDS020261328 



The He flux map is peaked and faintly extends toward 
the South. The velocity dispersion map has a peak around 
120 km 5"' in the South-East. The velocity field shows a clear 
but irregular velocity shear of about 100 km 5"' that extends 
over 0.8" in radius. Fitting a rotating disk model suggests that 
only the central part of the galaxy is detected since the plateau 
is not reached. It also shows that the velocity dispersion peak 
cannot be explained by the beam smearing effects. For the 
same reasons as VVDS020147106, VVDS020261328 has a 
dynamical mass (19.6 x lO'^'M©) 30 times larger than the stellar 
mass (0.6 x IO'^'Mq) suggesting that dark matter and cool gas 
dominate the gravitational potential of this galaxy. The SFR is 
the lowest of our sample (see Table |3]l. This galaxy is classified 
as a pertubed rotator mainly rotation supported (y„,flj^/cro = 3.5). 

VVDS220596913 

The Ha flux map shows two main components separated by 
12.5 kpc (1.5") that are not clearly distinguished in the /-band 
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Fig.A.5. Similar as Figure lATTI for VVDS220596913. The North feature is compatible with sky line residuals, /-band map is a 
CFH12K/CFHT image. 
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Fig. A.6. Similar as Figure Ia3] for VVDS 2205 84 167. The North- West feature is compatible with sky line residuals. 



image. The brigthest component is located in the western side 
and is composed by two peaks of equal intensity. The eastern 
component is much fainter. It has a smooth velocity gradient and 
two peaks in the velocity dispersion map (around 100 km s^^). 
One corresponds to the brightest Ha knot, and the other is more 
diffuse and corresponds to the transition region between the two 
components where the Ha emission is the faintest (signal to 
noise ratio around 4). The high velocity dispersion between the 
two components can be a signatur e for m erging as for Hickson 
compact group H31 ( Amram e t ai]|2007l) . The velocity shear of 
the faint component is misaUgned by around 20° with respect 
to the one of the brigth component suggesting, as does the flux 
distribution, that this system is composed of two galaxies in the 
process of merging. The velocity separation of the two com- 
ponents is estimated to ~ 150 km s"'. The merging timescale 
would be around 1.2 Gyr. Although a large and unique rotating 
disk c onsituted by large clumps as seen in iBournaud et al.l 
(|2008|) cannot be excluded (leading to Vmax/o-Q - 2.3), fitting 



the brightest component alone gives better results in terms 
of suggesting that the main component is rotating with a 
maximum velocity of around 200 km s^^. We fixed the center 
between the two flux peaks of the main component, and used 
an inclination of 55°. Assuming that both components have 
the same dimensions (half the whole system dimensions), we 
obtained a dynamical mass of around 8 + 2 x 10"^ (domi- 
nated by rotation) for the bright component and 4 + 2 x lO"* M© 
(dominated by dispersion) for the faint one. This imply a total 
dynamical mass of ~ 12x10'" Mq of the same order as the one 
computed considering a unique object (~ 13 x IO'^Mq). This 
object is classified as a major merger system. 

VVDS220584167 

VVDS220584167 is the object that presents the most 
extended Ha map of the sample. It has an elongated peak in the 
velocity dispersion map close to the center. The observation is 
affected by a strong sky line residual (16 195 A, see Figure |8]l 




Fig. A.7. Similar as Figure|A3]for VVDS220544103 
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Fig. A.8. Similar as Figure|A3]for VVDS220015726. 



that induces large uncertainties on line measurements, mainly 
on the southern side. The velocity field shows asymmetry and 
the galaxy presents a distorted morphology in the /-band image 
which makes possible the existence of a stro ng bar (see for 
comparison UGC 08937 in lEpinat et all 12009). Moreover, the 
observed displacement between /-band and Ha morphologies 
suggests a violent star formation event compatible with the high 
SFR of this galaxy. The velocity field is reasonnably fit by a 
simple rotational disk model, and the central velocity dispersion 
peak is not accounted for by the model since it does not simulate 
the effect of a bar. VVDS220584167 is the most massive disk 
of the sample with a dynamical mass of 21.1 x 10'"Mq. It 
is also the most rotation-dominated object {Mg/M,j- = 18 and 
Vmax/o'Q = 5.9). Thus this galaxy is classified as a rotating disk 
with a bar 

VVDS220544103 

The Ha flux map of VVDS220544103 shows two peaks 
separated by around 6 kpc (0.75"). It also shows a distorted 



morphology both in the line flux map and in the /-band image. 
A faint velocity gradient (~ 30 km i"') is observed toward 
the main northern peak, whereas a strong velocity gradient 
(~ 120 km s^^) is observed at the South of the fainter southern 
peak. The two gradients are in the same direction and the 
velocity field is continuous accross the whole object. The 
isovelocities contours are perpendicular to the distorted and 
elongated morphology. The strongest peaks in the velocity 
dispersion are located at the edges where the signal to noise 
ratio is low. Excluding these peaks, the velocity dispersion map 
shows an elongated peak close to the faint Ha flux peak. Fitting 
a simple rotating disk model to the whole system leads to a 
ratio Vmax/cro - 2.1. However, this object is hard to fit with 
such a model and the dynamical mass estimates just give an 
indication on its total mass. From our data, it is not possible to 
separate two components and compute mass estimates for both. 
We estimated the total dynamical mass to 10.3 x IO^'^Mq using 
a unique rotating disk object hypothesis. However, one might 
guess that at least one component could be seen neai^ly face-on 
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Fig. A.9. Similar as Figure|A3]for VVDS220014252. 



and thus have a larger maximum velocity and thus a larger total 
mass. The velocity separation between the two components is 
around 40 km '. The merging timescale for this system would 
be around 0.5 Gyr. This object could be two or more galaxies in 
a merging process and is classified as a merging system. 

VVDS220015726 

The central regions give a clean signal, while the low signal 
regions are affected by sky residuals due to the proximity of two 
strong sky lines next to Ha (15053Aand 15056A). The Ha map 
is peaked at the center and it has also a peak at the same location 
in the velocity dispersion map (~ 100 km 5"'). The kinematical 
maps of this galaxy are well reproduced by the rotating disk 
model. In particular, this model shows that the central velocity 
dispersion is due to beam smearing effects (see Figure [TTI). The 
relatively high;^^^ value of the fit is only due to the fact that this 
object shows the highest velocities. VVDS220015726 is clearly 
dominated by rotation iV,„ax/o'o = 8.5) since it is the fastest 
rotator of the sample (V„,ax = 323 km but also the object 
with the lowest velocity dispersion (ctq = 38 km i"'). Its total 
dynamical mass (16.7 x 10'*^ Mq) is compatible within the error 
bars to its stellar mass (6.2 x 10"^ M©). This indicates that the 
halo contribution is low in the central parts of this galaxy. This 
object is classified as a rotating disk. 



150 km s~^) are also observed at the North-East edge, with a 
rather good signal to noise ratio (larger than 5). The velocity dis- 
persion is not peaked at the center of the galaxy but is high ev- 
erywhere. Moreover, the mean velocity dispersion of this galaxy 
is the highest of the whole sample (92 km 5"'). The model shows 
that the velocity dispersion is not due to beam smearing effects 
(see Figure [TTI). The broad and asymmetric profiles suggests the 
possibility to h ave intrinsic double profiles that can be signatures 
of interactions (Amram et alj|2007 ). Moreover, the kinematics is 
not dominated by rotation since Vma.Jo'o = 1-1 and the total dy- 
namical mass (10.2 x 10'" M©) is dominated by the dispersion 
mass (7.9 x 10'" Mq). For these reasons, this galaxy showing 
the highest SFR of the sample is classified as a perturbed rotator 
possibly in the process of a minor merging event. 



VVDS220014252 

The Ha flux map of VVDS220014252 shows an elongated 
peak and diffuse emission in the outer parts, compatible with 
the /-band morphology. The peak does not perfectly match the 
center of external isophotes. The Ha flux map also suggests the 
presence of an arm in the western side. The velocity field is the 
one of a rotating disk except in the eastern side where it shows 
unexpected high velocities. A rotating disk model correctly fits 
the velocity field. The fit is better (smaller x^) when the cen- 
ter does not match the Ha flux peak, but the center of external 
isophotes. Adopting an inclination of 56°, the deprojected max- 
imum rotation velocity is quite low (103 km i"'). The model 
shows that the maximum rotation velocity is reached close to 
the kinematical center. The North- West side of the velocity field 
shows a bump and thus a small decrease. Around this bump, 
profiles are broad and asymmetric. Broad profiles (more than 



